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ABSTRACT 


This monograph presents a detailed analysis of the six- 
dimensional shapes of a particular type of Gravity Gradient 
boom bent and twisted by the thermal stresses induced by a 
solar thermal field. Contained within the analysis is a 
general method of solution, particularly suited for the de- 
flection analysis of very long thin-walled members of open 
section having a forcing function that varies along the length 
and depends on the position and orientation of the cross 
section relative to a fixed reference frame. 

The results given are intended to illustrate the effects 
of transverse -torsional coupling on the static -thermal- 
equilibrium shapes of the boom. It is shown that, for the 
Gravity Gradient boom studied, the thermal equilibrium 
shape is not unique and that it is not apparent which of the 
possible shapes are stable. It is thus evident that a planar 
assumption for thermal bending of any similar Gravity 
Gradient boom must be strongly justified before it can be 
assumed credible. 
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THERMAL BENDING PLUS TWIST OF A THIN-WALLED 
CYLINDER OF OPEN SECTION WITH APPLICATION 
TO GRAVITY GRADIENT BOOMS 

by 

Harold P. Frisch 
Goddard Space Flight Center 


INTRODUCTION 

With the advent of Gravity Gradient type satellites, the problem of predicting the motion 
of extremely long appendages has become important. These appendages are commonly referred 
to as "booms” and can usually be classified under the broad title of thin-walled cylinders of 
open cross section. 

Before any attempt can be made to accurately approximate the thermal bending of a boom and 
incorporate it into a dynamic model of a Gravity Gradient satellite one must first understand what 
the static thermal equilibrium shape of a boom bent and twisted by thermal stresses will be. 

The standard assumption made and used to obtain a first-order approximation of thermal 
bending is to treat the boom as a seamless cylinder of non-symmetric cross section bent by a 
definable thermal-stress distribution. The work here discussed treats the cylinder as one of open 
section; hence, it can take into account the effects of thermal torque and transverse-torsional 
coupling on thermal bending. It is shown that when the boom is treated as a cylinder of open 
section, more than one thermal- equilibrium shape may exist and that the trivial solution to the 
derived set of equations is the zero-twist solution; this is just what would be predicted from the 
seamless-cylinder assumption. 

In order to obtain a numerical solution to the derived equations, a worst-case estimate of the 
thermal stress distribution about the cross section of a silver-plated DeHavilland-type boom is 
derived. This is done by stating a set of reasonable assumptions and deriving from them an equa- 
tion that defines the temperature distribution for any cross-sectional orientation relative to the 
sun line. 

Since the three positional coordinates of the points on the longitudinal axis of the boom are not 
enough to define the cross-sectional orientation, three rotational coordinates must also be supplied; 
that is, an additional coordinate system traveling along the booms’ length with its axes aligned with 
the principle axes of inertia of the cross section. This system is related to the fixed system by a 
set of Euler angles, and the resulting equations of bending and twist are written in terms of them. 
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By writing the deflection equations as such, large-angle deflection relative to the fixed reference 
frame can be studied; because, relative to the local reference frame, the small-angle approxima- 
tion of bending and torsion has not been violated. 

The existence of a non-uniform temperature distribution implies that thermal stresses are set 
up which try to bend, twist, and elongate the boom. If we neglect the effect of elongation and make 
the Bernoulli-Euler assumption that plane sections remain plane in bending, we can derive the 
longitudinal stress distribution and write it as a function of the components of the resultant thermal- 
bending-moment vector. 

In general, the existence of longitudinal stresses implies the existence of shear stresses. It 
can be shown that they are mathematically related through a partial differential equation. The so- 
lution to this equation is shown and the resulting shear- stress distribution is derived. This is 
used to determine the resultant thermal torque that tends to twist the boom about the longitudinal 
axis. 

The method used to solve this particular problem of bending plus twist is derived by extending 
the results of S. Timoshenko (References 1 and 2) to take into account the constraints of this prob- 
lem. He showed that the problem of thin-walled members of open section bent and twisted by a 
given force distribution was separable. That is, if the deflection is small-angle, the equations of 
bending and torsion are independent and may be solved separately. 

The prime constraints which must be accounted for in the thermal bending problem are that a 
stress distribution rather than a force distribution is given and that the deflection will be large- 
angle relative to a fixed reference frame. These facts violate the initial assumption of Timoshenko T s 
problem; hence, his conclusions are not directly applicable. 

When a stress distribution is defined and large-angle deflection exists, the bending and torsional 
equations are coupled and must be solved simultaneously. 

It will be shown that the solution is given in terms of a set of simultaneous differential equations 
having all their boundary conditions but one defined at the root. This missing condition is defined 
at the tip and is dependent upon whether or not the tip is free to warp. 

In order to solve the set of equations numerically, the missing condition at the root must be 
guessed at and the resulting solution examined to see at what points the end condition is satisfied. 

At each point where the boundary condition is satisfied, a solution is said to exist and referred 
to as a "solution length." It was found that a particular assumed initial condition can satisfy the 
boundary condition at a number of points along its path of integration. Conversely a boom of a 
given length can have more than one initial condition that satisfies its boundary condition. This 
fact leads to the idea of more than one thermal -equilibrium shape. 

The results of this report show how the solution length changes as the missing initial condition 
is varied over its full range. By recording tip conditions at each solution length, it is shown how 
the tip deflection, tip twist, and direction of tip deflection change for each thermal- equilibrium 
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shape. The entire deflection pattern of a 100-ft boom is also shown for selected thermal equilibrium 
shapes and sun orientations. 

In order to interpret these results correctly it must be remembered that the bending and torsion 
equations are coupled and that the direction of deflection at any point depends on past history as 
well as local stress conditions. 

Furthermore, this report contains a general algorithm that may be used to solve similar prob- 
lems. The mode of presentation is designed to help the reader determine where he must make 
modifications in order to incorporate the constraints of his particular problem. 


SYMBOL LIST 


[A] 
BM X ( Z ) 
BM y ( z ) 
C 

c . 

E 

e 

e c 

G 

[T. T.E] 

[ij. jj. k,] 


Js 

K 


L 

n(s) 


coordinate transformation matrix 

thermal bending moment about X 2 ( z ) body axis at z 

thermal bending moment about Y 2 (z) body axis at z 

torsional rigidity 

warping rigidity 

Young T s modulus of elasticity 

distance between geometrical center and shear center of cross section 
coefficient of thermal expansion 
shear modulus 

orthonormal set of basis vectors parallel to axes of body triad [x 2 (z), Y 2 (z), 
z 2 ( z )], respectively, but having their origin at z = 0 

orthonormal set of basis vectors fixed in inertial triad and parallel to 
(x 1? Yj, zj respectively 

geometrical moment of inertia about X 2 (z) body axis 

geometrical moment of inertia about Y 2 ( z ) body axis 

solar radiation intensity 

thermal conductivity 

total length of boom (solution length) 

unit vector in direction of surface normal at position s 
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P = total perimeter of cross section 


Q(z), V(z) = quantities used to determine thermal torque coefficient 
r = radius of cross section 

R(z) = magnitude of deflection of point z on boom measured in X 19 Y x plane 

s = arc length measured from outer to inner seam around cross section 
SL = unit vector parallel to sun line 

t = thickness of cross section 

T( s , z ) = temperature (absolute) at point (s, z) on boom surface 

T m (z) = mean temperature of cross section at z 
T 0 = ambient temperature 

T c (z) = thermal torque due to shear stress distribution at z 

Xj = axis parallel to principal axis of inertia at z = 0 oriented so that (x x , y 19 Z 1 
is right-handed system, outer seam in positive X x region for X = -1, outer 
seam in negative X 1 region for k = +1 

(X x , Y x , z x ) = inertial triad 

X 2 ( Z ) = principal axis of inertia normal to Y 2 (z) in plane of cross section, outer 
seam in positive X 2 ( z ) region for k = -1, outer seam in negative X 2 ( z ) 
region for \ = +1 

x 2 (s) = X 2 (z) coordinate of point s on cross section 

[x 2 (z), Y 2 (z), Z 2 (z)] = body triad at z, principal axes of inertia 

Yj = axis parallel to mass symmetry axis of cross section at z = 0, k = +1 or 
-1, overlap in negative Y x region 

Y 2 ( z ) = mass symmetry axis of cross section, overlap in negative y 2 ( z ) region 
for k = +1 or -1 

y 2 (s) = Y 2 (z) coordinate of point s on cross section 

Z x = axis tangent to centroidal axis at z = 0, positive in direction of increasing z 

Z 2 ( Z ) = principal axis of inertia tangent to centroidal axis at z in direction of in- 
creasing z 

z = arc length measured along centroidal axis from root to tip 
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absorptivity 
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e 

e z (s, z) 
0 

6>(z) 

[ 0 1 ( Z )’ ^2 < Z )’ < P ( Z )] 
A. 

X 

7T 

-j - a(z) 
P( z ) 

Px O)' Py ( z )' Pz ( z ) 

O* 

°z (s. z) 

r( s , z ) 

res) 

4> 

cp(z) 

¥ 

'/'o’ a o 
«*) 


= distance between geometrical center and centroid of cross section 

emmissivity of surface 
thermal strain at (s, z) 

= direction of bending if zero twist bending existed, normal to neutral axis 

direction of deflection, i.e., inclination of R(z) from negative X x axis meas- 
ured in a right-handed sense 

Euler- angle rotation sequence used to define orientation of body triad 
relative to inertial triad 

+1 if direction of increasing s is clockwise about positive Z 1 axis in right- 
handed sense 

-1 if direction of increasing s is counterclockwise about positive Z 1 axis in 
right-handed sense 

angle between sun line and Z 2 (z) axis at z 
curvature vector at z 

components of curvature vector at z relative to body triad 

radiation coefficient 

longitudinal thermal stress at (s, z) 

shear stress at (s, z) 

fl for SL * n(s) ^ 0 
1.0 for SL * n(s) < 0 

angular amount of overlap 
boom twist at z 

direction of bending if "in plane" bending existed 
values of 0(z) and a(z) respectively at z = 0 

angle between component of sun line in cross-sectional plane at z and the 
negative Y 2 (z) axis, measured in direction of increasing s 
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TEMPERATURE DISTRIBUTION 


The foundation of any analysis of thermal bending is the accurate description of the tempera- 
ture distribution in the specimen under study. An extremely accurate description, however, in- 
evitably leads to equations that can be solved only when simplified. 

The thermal equations presented in this section attempt to give a worst-case estimate of the 
thermal gradients; hence, the second-order effects that tend to reduce the thermal gradients are 
neglected. By considering the worst-case estimate, this report will attempt to define the maximum 
limits of deflection and accent any effects that would not be predictable from a less detailed analysis. 

The thermal equations used and stated below were derived by simply equating heat-in to heat- 
out for an element on the outer face of the cylinder 0 £s £ 2 ttt and for an element on the overlapped 
portion of the cylinder 2ttt <s £ P, where: 

s = arc length measured positive in a circumferential direction from the outer seam 

r = cylinder radius 

P = total perimeter of the cross section. 

These equations are based on certain assumptions: 

1. Heat is conducted in a circumferential direction around the entire perimeter 0 < s < p . 

2. The entire cross section 0 £s < P loses heat by radiation. 

3. The radiant heat absorbed by an element on the sunlit side of the cylinder is proportional 
to the cosine of the angle between the surface normal at the element and the sun’s rays. 

4. The radiant heat absorbed by an element in the overlapped portion of the cylinder is pro- 
portional to the heat radiated by the element directly above it on the outer face. 

5. The radiant heat absorbed by an element due to internal radiation is assumed zero, since 
its inclusion would reduce the thermal gradients. 

6. There is zero contact between the overlap and the overlapped portions of the cylinder. 

Hence, heat cannot be conducted by point contact, and the effect of friction and stiction may 
be neglected in the analysis of bending and twist. 

By directly comparing the assumptions stated with the following equations one can be assured 
of their validity 


Kt 


d 2 T(s ) 
ds 2 


o-e T 4 (s) - J s a s T(s) [SL • n(s)] 


for 01s < 2ttt and 


Kt 


d 2 T(s) 
ds 2 


creT 4 (s) - a s T 4 ( s - 277T ) 
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for 2ttt < s < P, where 


dT(s) 

ds 


s-0 

s=p 


0 


The following symbols have been used in the above equations: 
t = thickness of cross section 
K - thermal conductivity 

T( s ) = absolute temperature at position s 
a - radiation coefficient 
e = emissivity of material 
J s =- radiation intensity of the sun 
a s = absorbtivity of material 
SL = unit vector in direction of sun line 

n(s) = unit vector in direction of surface normal at position s 
fl for SL • n(s) > 0 

T{ s) - 1 _ 

lO for SL n(s) < 0. 

It will be noted that the solution to the above equation is a boundary value problem. However, 
it is of a type which can be easily solved numerically with great accuracy by using a steepest- 
descent technique. 

COORDINATE SYSTEM 

For a thin-walled cylinder of open section (a boom) clamped at the root, the cross section may 
be oriented in one of two different ways (see Figure 1): 

1. Such that the path from the outer seam to the inner seam is traced by a counter clockwise 
rotation of the radius vector or, 

2. Such that the path from the outer seam to the inner seam is traced by a clockwise rotation 
of the radius vector. 
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Y 1 Y, 

CASE I CASE H 

COUNTERCLOCKWISE SEAM ROTATION CLOCKWISE SEAM ROTATION 


Figure 1 —Inertial axes orientation. 

The two coordinates (s, z) will locate any point on the surface of the boom. The coordinate z 
measures arc length along the centroidal axis of the boom measured positive from the clamped root. 
The coordinate s measures arc length around the cross-sectional perimeter, at the point z , measure 
positive from the outer seam. At the inner seam for both cases s = P,the cross-sectional perimeter 

The constraints of the thermal-bending problem require that the cross-sectional orientation at 
every point be known and that large-angle deflection relative to a fixed reference frame be predict- 
able. Since the equations of bending and twist relative to a local reference frame do not violate the 
small-angle assumption of deflection, an additional coordinate system is defined which travels 
along the boom’s length. 

At any point z along the boom’s length, a coordinate system [x 2 (z), Y 2 ( z ), z 2 (z)] can be de- 
fined. These three axes form a right-handed orthogonal coordinate system parallel to the three 
principal axes of inertia with the origin along the centroidal axis. As the coordinate z increases 
from 0 to L (the total boom length) and the boom is deflected, the triad [x 2 (z),y 2 (z),Z 2 (z)] 
will be translated and rotated with respect to an inertial triad [x t , Y lf Z t ] . The inertial triad is 
defined to be coincident with the body triad at z = 0, the root. 

At any point along the boom’s length z, the axes of the body triad [x 2 (z), Y 2 ( z ), z 2 ( z )] are 
defined as follows (see Figure 2): 

Y 2 ( z ) = Mass symmetry axis of the cross section directed downward so that the overlap is in 
the negative y 2 ( z ) region for both cases. 

X 2 (z) = Principal axis of inertia in cross-sectional plane normal to Y 2 (z) and through centroid. 
In case 1 (Figure 1) for counterclockwise seam orientation, the outer seam is in the 
positive x 2 (z) region. In case 2 the outer seam is in the negative x 2 (z) region. 
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Z 2 (z) = Axis parallel to the centroidal 
axis at z, with origin at the cen- 
troid of the cross section and 
directed in the direction of in- 
creasing z. 

The orientation of these three axes rela- 
tive to the inertial triad, at any point z along 
the boom T s length, can be completely defined 
by three successive Euler-angle rotations. 

Let i 1? T 1? kj be a basis of three ortho- 
normal vectors fixed in the inertial triad and 
parallel with the X l9 Y l9 and Z x axes, respec- 
tively, with their origin at the centroid of the 
cross section at the root z = 0. 

Let i,J, k be another basis of three ortho- 
normal vectors parallel to the axes of the body 
triad X 2 (z) , Y 2 (z), and Z 2 (z), respectively, 
but having their origin at the centroid of the 
cross section at the root z = 0. 

By defining the angles & 1 (z), & 2 (z), and 
<p(z) as the three successive Euler-angle rota- 
tions which will define the orientation of the 
body triad relative to the inertial triad at the 
point z, as shown in Figure 3. It follows that 
the vectors , Tj , k 1 can be written in terms 
of the vectors I, j*, k by means of the trans- 
formation equation 




i 


i 


1 0 

0 cos(9 1 (z) 

0 sin^ (z) 


0 

-sin 6^ (z) 
cos& 1 (z) 


cos 0 2 (z) 
0 

-sin &2 (z) 


0 

1 

0 


sin 6 2 (z) 
0 

cos 6 2 (z) 


cos cp(z) 
sin cp(z) 
0 


-sincp(z) 
cos qp(z) 
0 



sintfj s 
-cos (9j ; 


-cos$ 2 sin 9 

-sin 0 1 sin d 2 sin qp + cos & 1 cos <p 
cos Q x sin 6 2 sincp + sin 6^ cos qp 


sin # 2 

-sin 6 x cos 0 2 
cos 6 X cos 6 2 



i 


[A] 


i 

k 
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Since [A] is an orthonormal transformation matrix, its inverse is equal to its transpose and 
hence 



At any point z along the boom's length the unit vector k* tangent to the centroidal axis is de- 
fined by the vector equation 


_ dX 1 (z) ^ dY j ( z ) dZj ( z ) 

k* dz *1 + dz -j l + dz ^1 * 

where [x t (z), Y x (z), Z x (z)] are the inertial coordinates of the point z on the centroidal axis of 
the boom. By definition the vector k is parallel to k*; hence, from the transformation equation it 
follows that 

k* = k r sin# 2 (z)i 1 -sin# 1 (z) cos 6 2 (z) j x + cos 0 x (z ) cos 0 2 (z) kj . 

Set the components of the vectors k* and k equal to each other and solve the resulting differ- 
ential equations; the inertial coordinates of any point z along the centroidal axis can then be obtained; 
that is: 


dX, (z) 


ii 

N 

sin<9 2 (z) 

dYi (z) 


dz 

- sin^j (z) cos 9 2 (z) 

dz, (z) 


dz 

cos 9 j (z) cos 9 2 (z) , 


where 


dXj (z) 

i 

/*— N 

N 

' 

1 

ft 

N 

dz 

'z = 0 ' dz 1 

z = 0 dz 1 


The differential equations defining 6 X (z) and 0 2 (z) will be derived in the section entitled "Equa- 
tion of Bending." 
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THERMAL STRESS 


The formulation of the equation defining thermal stress will be based on the Bernoulli -Euler 
assumption that plane sections remain plane in bending and that the effects of lateral contraction 
may be neglected. 

It is shown in the theory of thermo elasticity that Hooke’s Law must be extended to take into 
account the effects of thermal expansion (see Reference 3). Thus, the thermoelastic equation for 
longitudinal thermal stress at a point on the surface of the boom is given by 

cr z (s, z) = E (s, z) -e c [T(s, z) -T 0 ]j ■ 


where 


(S' z ) 
(s, z) 

e c 

T(s, z) 

T 
1 o 

E 


longitudinal thermal stress 
longitudinal thermal strain 
thermal expansion coefficient 
absolute temperature 
absolute ambient temperature 
Young’s Modulus of Elasticity. 


It should be noted that there is no restriction made on the way the temperature T(s, z) varies 
around or along the boom. The Bernoulli -Euler assumption that plane sections remain plane in 
bending, however, requires that the displacement u(s, z) of any point on the cross section be given 
by a linear function of its x 2 (s), y 2 (s) body coordinates. Hence, the longitudinal strain must be 
of the form 


du(s y z') 

e z ( S ’ z ) “ ““Vz “ f 0 ( Z ) + X 2 O) f l ( Z ) + y 2 ( S ) f 2 ( Z ) * 

where f 0 (z), (z) and f 2 (z) are obtainable from the equilibrium conditions. These equilibrium 

conditions require that the resultant force and moments about the principal axes of the cross 
section vanish, that is 


f 

Jo 


ta z ( s , z) ds 


■f 

•'o 


to-, (s, z) y 2 (s) ds 


f 

Jo 


t<? z (s, z) x 2 (s) ds 


If we directly substitute e z (s, z), as defined above, in the longitudinal- thermal- stress equation 
and make use of the equilibrium conditions, f 0 (z), fj(z), and f 2 (z) can be derived. It 
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follows that 


BM y (z) BM x (z) 

(s, z) = Ee c [T b (z)-T(s, z)] - x 2 (s) + y 2 (s) 


where 


f 

J o 


Iv - I t y 2 2 (s) ds 


f 

J o 


«- Y . t x 2 2 (s)ds 


T -(o 


- \ f T ( s ’ z ) ds - TF f 
•'0 *'0 


T(s, z) ds 


•'o 


BM x (z) - e c E t T(s, z) y 2 (s) ds 


BM y (z) = ‘ e c Et 


JV 

•'0 


z) X 2 (s) ds 


SHEAR STRESS 

It is shown by Timoshenko in Reference 2 and in most books on elasticity that for thin-walled 
members the shear stress r( s, z) is related to the longitudinal stress < j z (s, z) by the partial dif- 
ferential equation 


<JtT(s, z) aa 2 (s, z) 

— h — = " t — Tz — • 

By directly substituting the equation derived for thermal stress in this equation, it follows that 


(9tr(s, z) 
d s 


tEe 


~ aT m ( Z ) _ 3T(s, z)’ 
dz dz 


tx 2 (s) 5BM y (z) 
I <9z 


ty 2 (s) <9BM x (z) 
I dz 


The elements of this equation evaluated at point z depend on the temperature distribution, 
which depends on the relative sun position; this, in its turn, depends on Euler angles 6 1 (z), d 2 (z), 
and cp(z). Therefore, in solving this problem numerically, it is convenient to differentiate with 
respect to the angle of twist cp(z) rather then z . 
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If we apply the chain rule of differential calculus and integrate both sides of the equation from 
0 to s, the expression defining shear stress reduces to 


tr(s, z) 


d(p(z) 
* dz 


f S /„ p T m ( z > dT(s, z)] * 2 ( s ) dm y ( z ) y 2 < ! 

Jo \ Eec <3<p(z) <M z ) I y acp(z) + I X 


s) <9BM X (z) 
dy(z)~ 


ds , 


where the functions x 2 (s) and y 2 (s) that define the coordinates of the point s on the cross section, 
in the x 2 (z), Y 2 (z), Z 2 (z) coordinate system, are different for the case of clockwise rotation and 
for the case of counterclockwise rotation of 
the cross section, as seen from Figure 4. 

From the geometrical relations shown in this 
figure it can be seen that 


X 2 (s) 

, . /s 

4 >\ 

A. r sm l — 

" l) 

y 2 ( s ) 

/s 

<t>\ 

- r cos l — 

" 1) 


+ s , 



I x 2 (») 


CASE I 


CASE H 


where Figure 4— Body axes orientation. 

8 = distance between centroid and geometrical center of cross section 

. 4> 

r sin 2 

~ 4> 


77+^ 


-g: 


if clockwise seam rotation, Case n 
if counterclockwise seam rotation, Case I. 


Furthermore, from the equations defining the geometrical moments of inertia it follows that 

I x « t r 3 £(77 + + sin - (277 + 0)tr8 2 

r 3 ^7T + |) - sin y] • 


= ‘ 


Substituting the equations for x 2 (s) and y 2 (s) in the expression for shear stress and performing 
the integration indicated gives 


tr(s, z) 


- t 


dcp( z 
dz 


) | f ( z ) f dT(s, z) *1 Xr 2 (z) |" 0*] 

-S Ee c [ s wry - J o a«K«) ds J- t; [~ cos v 7 " ^ + cos 2\ 


r 2 <9BM x (z) 

1. a<P(z) 


.(?-« 


+ • ^ 

+ sin o ~ 


. <£ 

s sin ~2 

yTf) 
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This expression defines the shear stress 
at any point on the boom’s surface. It is a 
vector having the units of force and directed 
so that it is normal to the longitudinal axis at 
z and tangent to the surface at (s, z) in the di- 
rection of increasing s, see Figure 5. 

BENDING PLUS TWIST 

For thin-walled members of open cross 
section there is an axis along which any 
applied transverse force distribution will pro- 
duce pure bending with zero twist. This axis i 
its intersection with a cross-sectional plane as 



COUN TERCLOC KWI SE C LOC KWI SE 

SEAM ROTATION SEAM ROTATION 

CASE I CASE H 

Figure 5— Direction of shear-stress vector. 

; referred to as the "shear center axis" and 
the "shear center." 


For the case of small-angle deflection relative to a fixed reference frame, the deflected shape 
of a member of open cross section when subjected to a force distribution along any axis can be 
divided into two independent parts and the results combined, that is: 


1. The determination of bending for that force distribution applied along the shear center axis, 
and 

2. The determination of twist resulting from the moment of that force distribution about the 
shear center axis. 


As mentioned in the introduction, this method has been discussed and proved by S. Timoshenko 
in Reference 2. 


The solution to the thermal-bending problem involves taking into account various constraints 
not considered by Timoshenko in the solution to the above problem. These are: 

1. Relative to a fixed reference frame, the bending and twist of the cylinder of open section 
cannot be assumed to be small-angle over the entire length; hence, the standard small-angle 
assumptions made in his analysis are not valid for the thermal-bending problem relative to 
a fixed reference frame. 

2. The thermal stress distribution at any cross section is a function of the sun line orientation 
relative to the principal axes of inertia of the cross section. Hence, three additional equa- 
tions are needed to define cross-section orientation relative to the fixed reference frame. 

3. When the cylinder is twisted, the principal axes of inertia at the root are not parallel to the 
principal axes of inertia of a cross section displaced from the root; hence, the bending 
stiffness at the displaced cross section measured about the axes fixed at the root is not 
equal to the bending stiffness at the root measured about these same axes. 
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In order to solve the thermal-bending problem, a local coordinate system has been defined at 
every point along the centroidal axes. Since curvature at every point can be measured with re- 
spect to these local axes, the standard small-angle assumption of bending and torsion may be em- 
ployed; Timoshenko’s results are applicable over every element of length, relative to the local 
axes. Thus, the equations of bending and torsion are independent over each element of length and 
may be solved separately. The orientation and position of the local axes relative to the fixed axes 
at each integration step along the boom’s length may be determined by employing the derived co- 
ordinate transformation equations. 

Besides defining the boom’s shape in a fixed coordinate system, this information is needed to 
calculate the appropriate thermal stress distribution, which is a function of the relative sun posi- 
tion at each integration step. 

THERMAL TORQUE 

S. Timoshenko in Reference 2 shows that for a thin-walled cylinder of open cross section, an 
applied torque will be balanced partly by the cylinder’s torsional rigidity and partly by the cylinder’s 
warping rigidity. This is expressed by the following formula: 

T sc ( z > = " C 1 p"' ( z ) + c< p' O) , 


where 

T sc (z) = applied torque 
C x - warping rigidity 
C = torsional rigidity 

<p(z) = twist, positive about positive z 2 (z) axis. 

The boundary conditions associated with these equations are 

<P( Z )I z =0 = 0 . 

<P* <«)l,-o = 0 , 

and 


tp" ( z )L=l 


0 , if the tip is free to warp, 


<P* (Z)I Z=L 


0 , if the tip is not free to warp. 
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It should be noted that the rate of twist about the Z 2 (z) body axis as derived from the Euler 
rate equations is cp' (z) + (z) sin Q 2 ( z ) . However, the factor 6^ ( z ) sin0 2 (z) is assumed to be 

approximately equal to zero; hence, the rate about the Z 2 (z) body axis is accurately given by <p' (z). 
The assumption that this factor is approximately zero follows from the Bernoulli- Euler assumption 
of bending, which states that the transverse bending moment vector lies in the plane of the cross 
section (see Section "Equation of Bending"). 

In Reference 2 it is further shown that the torsional rigidity C is given by 

C = 3- Pt 3 G , 

where 

G = modulus of elasticity in shear 
t = thickness of cross section 
P = total perimeter of cross section 

and the warping rigidity C x is given by 

J\ <t> ( . *\ <t>V 

2 o sm ~2 - \jt + ~2 J cos ^ 

? 

7T + ~2 ~ sincfc 

where 

<t> = total angular amount of overlap 
P 

r = 2 n +0 = radius of the cylinder of open section. 

The applied torque T sc (z) may be calculated by integrating the torque due to the shear-stress 
distribution around any axis normal to the cross section. It will be seen that if the shear center 
axis is chosen as the axis about which the torque is evaluated, the resulting expression reduces to 
an equation that is independent of the bending- moment components. 

For a cylinder of open section whose diameter remains constant along the length (as has been 
assumed), the distance e between the shear center and the geocenter of the cross section is given by 



This equation is also derived in Reference 2 and many other books on strength of materials. 
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The element of torque dT gc ( z ) due to the shear stress at an element ds of unit longitudinal 
length on the surface is 

dT sc (z) = xtV(s, z ) ds , 

where 

K = ± 1 as previously defined 

e is the vector from the shear center to the surface element ds . 

s 

The total torque T sc ( z ) is given by the integral of the above expression over the entire perimeter, 
i.e., 

T sc (z) = A. e g x tr(s, z) ds . 

•'o 

These vectors are shown in Figure 5 for both clockwise and counterclockwise seam rotation. 

It is evident from the figures shown that for both cases 

e s xtT(s, z) - e x tr( s , z) + r xtr(s, z) , 

where 

e is the vector from the shear center to the geocenter of the cross section 
r is the radius vector. 

From the definition of cross-product it follows that 

e x "tr( s , z ) = -|e tr(s , z ) sin ^ + ^ ("r “ k , 

which by the fundamental trigeometric identities reduces to 


e x tr( s , z) 


{• 


e tr(s, z) cos 



for both cases. 
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Dropping the vector notation used above (which is no longer needed since T sc is parallel to k) 
and directly substituting in the total torque integral gives 


T sc (z) - X. e j tr(s, z)cos 


I tr(s, z)cos (-p - ds + kr j tr(s, z)ds , 

J 0 


where 


tr(s, z) 


d(p(z) 
* dz 


Ee 


' dT m ( z ) 
s dcp(z 


z) r s 

> 'Jo 


<?T(s, z) 

dy(z) 


ds 


Xr 2 

«9BM y (z) 

’ ( s _ ±\ 

4>- 

T“ 

y 

0<P(z) 

- cos l r 2 j 

+ cos ~2 


r 2 (9BM x (z) 

— 



■ ^11 

sin 

/s <M 

. 0 

1 + sin -j ~ 

s sin 2 1 

I x d<p(z) 

\r " 2) 

r ( 77+ l)Jj 


In order to evaluate the above integral, the following definite integrals must be used: 



along with the shear center equation 


. <t> 

$ Sln 2 
cos 2 W 

7T + 2 

7 4> r . 

77 + ~2 ^ sin <p 
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Thus, after the appropriate substitutions and cancellations are made, the expression for total 
torque T sc (z) reduces to 


T 


sc 


<*) 


- \ tEe c 


d<P(z) 

dz 


„r . $ pi < 9T ».( z) 

r P [e s in 2 + 2\ dcp(z) 



r + e cos 



d l '(£, z) 

<?tp(z) 


d£ ds 


Let 


Q(s, z) 



T(<f, z) d# 


and note that 


T m ( z ) 


If 

•'0 


T(s, z)ds 


Then it follows that 


T S c ( z ) - - Ee c r J 0 {[esin 4 + 1] 5 J?(z) Z) " P [f cos (f - I) + i] 5 a%f^} ds 

This equation may be further condensed by making the substitution 

V(s,z) = T(s,z)ds-pJ* cos + lj Q(s, z) ds 

in the above expression. The total applied torque expression thus reduces to 


T sc O ) “ '^ Ee , 


dcp(z) <9 V(s, z) 
r dz 


acp(z) 


and the torque equation becomes 


or 


- UEe t r < * > ' ( z ) “ " Cj cp'" (z) + Cep' (z) 

t 

... r, t Ee c r av(p, Z )i , 

v (z) - |^k 2 + X. — dcp(z) J<P ( z ) - 0 , 
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where 


The method used in the determination of dV(P, z)/d<p(z) is not immediately obvious and will be ex- 
plained in detail in the next section. 

NUMERICAL EVALUATION OF BENDING MOMENT AND TORSIONAL TORQUE 

When solving the bending and torsional equations it is essential to determine how the tempera- 
ture distribution and hence the bending moment and torsional torque coefficients change with a 
change in relative sun position. Since it is computationally impracticable to do this at each inte- 
gration step in the solution, an alternate technique must be developed. 

Using the fact that the temperature distribution does not change radically for small changes in 
sun orientation and inclination, it is possible to compute the bending and torsional coefficients 
desired for any sun position by four-point linear interpolation when the bending and torsional coef- 
ficients are computed for a finite number of sun orientations and inclinations. 

By the appropriate definition of sun orientation V* and sun inclination a*, rectangular matrices 
of numbers defining the desired coefficients evaluated at incremental steps in y*, where 

0 < *p* S 2n , 


and at incremental steps in a*, where 


0 < a* < 

can be used for both clockwise and counterclockwise seam orientation if the appropriate sign 
changes are made. 

For both cases let (n/ 2) - a* be the angle between the positive Z 2 (z) body axes and the sun 
line. Furthermore, let V* be the angle between the negative y 2 (z) body axes and the component of 
the sun line in the cross-sectional plane measured positive in the direction of increasing s. Thus, 
for a particular sun position defined by [¥*, a*] , the identical temperature distribution will result 
for both cases when angles are measured as defined. 

It will be recalled that the temperature distribution was assumed to be given approximately by 
the solution to the following differential equation: 

Kt d ^ = ere T 4 (s, z) - J s a s 7(s)[sL * n(s)] 
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for 


0 < s <. 2jtt 


and 


\ 


d 2 T(s, z) _ 

Kt CT£ 


|t 4 (s, z)-a s T 4 (s - 2 t7t , z)j 


for 2ttt < s < P, where 


dT( s , 
3s 


z ) 


s-p 


0 


It follows from the definition of a* and w* just made that the unit vector SL in the direction of 
the sun line may be written in body coordinates as 

SL - K cos a* s in i - cos a* cos j + s in a* k 

(see under "Relative Sun Position"). 

In order to solve the problem proposed, the bending and torsional moments must be obtainable 
for any sun position relative to the body axes at any cross section. 

The following quantities are needed for the solution of the bending-plus -torsion equation and 
depend on the relative sun position. 


T m (z> 


■ if 
J 0 


T( s , z ) ds 


f 
J 0 


BM x (z) - e c E t T(s, z) y 2 (s) ds 


Q(s, z) 


f 

J o 


BM y (z) - - e c Et j T(s, z) x 2 (s) ds 


■if 


T(f, z)d# 


V(P, z) - |^e sin ^ T(s, z) ds - P J" cos (y - + lj Q(s, z) ds 
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The numerical value of each can be ascertained by solving the differential equations defining 
them simultaneously with the temperature -distribution equation and evaluating each at s = P, i.e., 


where 


dT „<s’ 

Z ) 

ds 


dBM x (s, 

Z ) 

ds 


dBM Y (s, 

Z ) 

ds 


dQ(s, 

_ Z ) = 

ds 


dV(s, 

_*) = 

ds 


z )| 

1 s=0 

= 


P T(s, z) 

e c Et T(s , z) y 2 (s) 

- e c E t T(s , z) x 2 (s) 
= p T(s, z) 


[e sin ^ + ^]t(s, z) - P jp cos (p - ^) + lj Q(s, z) 


T m (s, z)| s=o = BM x (s, z)| s=(j - B“y ( s -*)| bB0 = Q ( s ’ z )l s=0 = V ( S ’ Z )l s=0 


0 . 


Incrementally varying y* and a*, in an appropriate manner, gives a sufficiently dense distri- 
bution for each of the quantities, so that their values at any sun position can be accurately evaluated 
by interpolation. 

It was shown that the torsional torque was dependent upon the partial derivative of V(P, z) with 
respect to twist q>(z). In order to derive a matrix of numbers which defines this quantity, V(P, z) is 
numerically differentiated with respect to y* and the result is stored. This quantity is related to 
the desired quantity by the equation 

aV(P, z) aV(P, z) ay* 
a^p(z) ay* a<p(z) 

f 

A function 'P(z) analogous to ?* will be derived under "Relative Sun Position." From this expres- 
sion it can be shown that for small-angle deflection 

P(z) = 4> 0 - \<p(z) , 

where >A 0 is the orientation of the sun at z = 0. It is a constant defined initially; hence, 

dHz) x 
aqp(z) 
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But since '/'(z) is analogous to W* one may write 


<9V(P, z) ^ v <9V(P, z) 

<9cp(z) " ‘ K dV* 

If large-angle deflection is anticipated, the function ^(z) can easily be differentiated with respect 
to qp( z ) and must be introduced into the torsional equation. This must be done if booms several 
hundred feet long are to be studied. 

EQUATION OF BENDING 

Under "Coordinate System" it was shown that the Euler angles 6 > 1 ( z ) , d 2 ( z ) and cp(z) com- 
pletely define the orientation of the cross section located at distance z from the origin (z being 
measured along the boom T s longitudinal axis). 

From the elementary theory of bending it is known that the bending moment at a point along the 
boom T s length is proportional to the boom T s curvature at that point. Let ~p{ z) be the curvature 
vector at any point z along the boom T s length. In terms of the inertial triad, its magnitude and 
direction will be given by „ 


60 x ( z ) d$ 2 (z) d 0 2 ^ 

P( z ) = ” ' d^~ *1 + ~dz cos 6 y (z) jj + -g^ sin 6 X (z)k 1 . 


Making use of the coordinate transformation defined previously, it is possible to write this vector 
in terms of the body axes [T, k]; that is 


P(Z) - P x (z) i + P Y (z) j + P z (z)k . 

Since by the Bernoulli- Euler assumption the bending moment vector lies in the plane of the cross 
section, the k component of curvature p z (z) is approximately equal to zero; i.e., 


d6?i (z) 

Pz ( z ) “ di — sin # 2 (z) - 0 . 


(Note: this is the identical assumption made under "Thermal Torque.") 
Hence, 


P( z ) “ Px ( z ) 1 + Py ( z )T 
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where 


d 6 X (z) 69 2 (z) 

Px ( z ) = dz~ cos e 2 (z) cos <p(z) + sincp(z) 

d0j (z) d0 2 (z) 

Py (*) - dT~ cos 6 2 ( z ) sin<P(z) + — cos<p(z) . 

Thus 0 ( z ) and d 2 (z) can be obtained from the solution to the differential equations 

d0 i ( z > _ l r i 

“ di — - cos $2 (z) L -J °x < z ) cos <P(z) + p y (z) sin (p(z)J 

d0 2 (z) 

~~Si " Px ( z ) sintp(z) + p Y (z)cos cp(z) , 

where expressions for p x (z) and p Y (z) will be derived. 

At a particular cross section the thermal stress was shown to be 

* 

°- z (s, z) = E|e z (s, z)- e c [T(s, z)-T 0 ]| 

BM y (z) BM x (z) 

" Ee c [ T m ( Z )“ T ( S ’ Z )] " I y X 2 ( S ) + I x ^2 ( S ) * 

Hence the thermal strain, i.e., the longitudinal extension of the point x 2 (s), y 2 (s) on the cross 
section, is 


BM y ( z ) BM x (z) 

e z (s, z) - e c [T m (z)-T 0 ] - IyE x 2 (s) + — - E “ y 2 (s) , 

where the quantity e c [T m (z) -T 0 ] is analogous to the tension term that arises in the derivation of 
the equation of bending when bending and tension are considered (Reference 4, page 162). 

In order to determine magnitude and direction of bending, the effect of elongation will be 
neglected; i.e., it will be assumed that 


e ck( z )- T J - o • 

Hence the bending stress a* ( s, z) is given by 

BM y (z) bm x (z) 

a * ( s ’ z ) = IT - x 2 < s ) + y 2 ( s ) - 
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where positive bending moment components BM Y (z) and BM X (z) will cause positive bending in a 
right-handed sense about the Y 2 (z) and X 2 (z) body axes, respectively. 

Since bending moment is directly proportional to curvature, 

BM y (z) 

Ely = ^YO) 


BM X (z) 

EI X = " Px ( z ) - 

and the differential equations defining d 1 (z) and 6 2 ( z ) reduce to: 

6B j ( z ) 1 BM X ( z ) BM y ( z ) 

Si = cosS 2 (z) ,E I x cos <p(z) + EIy sintp(z) 


d& 2 (z) BM X (z) BM y (z) 

— ^ - - E I x sin <p( z ) + EI y coscp(z) . 

The bending-moment components as shown are functions of longitudinal arc length z . This is 
obvious since the sun position relative to the body axes changes as the boom is bent and twisted. 
Since the sun can be located by the angles ^(z) and a(z), the bending- moment components become 
functions of these two quantities, thus: 

BM x (z) = BM x (0,a) 

BM Y (z) = BM y (^ t a) , 

where 

~ - a (z) = the angle between the sun line and the longitudinal axis at z 

</f( z ) = the angle between the component of the sun line in the X 2 (z) , Y 2 (z) plane and the 
negative Y 2 (z) axis, measured in the direction of increasing s. 

These angles are easily obtainable by vector techniques and can be defined in terms of Q ± (z), 
e 2 ( z ), andcp(z), as will be shown. 
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RELATIVE SUN POSITION 


Under "Coordinate System” it was shown that the inertial triad [1 1 , ~j 19 kj could be written in 
terms of the body triad (i, j*, k) by the transformation equation 



and the unit vector k in the direction of the longitudinal axis is defined by 


k - sin $ 2 ( z ) i j - sin (z) cos 0 2 (z) j x + cos 6 x (z ) cos 0 2 (z ) kj 

written with respect to the inertial triad. 

Let 

SL = unit vector directed toward the sun 

< p 0 = ^( 2 >j 

'z=0 

a 0 = a < z >U- 

From Figure 6, showing the direction of the sun line vector SL, it can be seen that 

SL = cos a Q sin 4 s q i x ~ cos a Q cos j x + sin a Q k x . 

Since is measured in the direction shown in Figure 6 for the clockwise seam orientation case, 
and in the opposite direction for the counterclockwise seam orientation case, it follows that 

SL = k cos a Q sin \p 0 i 1 - cos a Q cos V^ 0 j x + sin a Q k x , 

where, as before, 

f+1 clockwise seam orientation 

\ = i 

L-l coimter clockwise seam orientation. 

From the definition of angle a(z) and the scalar product, 

~2 "* a (z) - cos' 1 [k * Sl] 


X, 



Figure 6 — Inertial triad showing the 
direction of the sun-line vector. 
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or 


a( z ) = ~2 “ cos” 1 [\ cos a Q s in i/> 0 s in 0 2 ( z ) + cos a Q cos i// () sin^ 1 (z) cos 6 2 (z) + sin a Q cos Q x ( z ) cos 6 2 ( z )] * 

In order to determine <A( Z ), write the sun vector in terms of the body triad, that is 

SL = KS X i — S 2 j + S 3 k , 


where 

Sj = cos a Q sin <A 0 cos Q 2 (z) cos cp(z) - K cos a Q cos 0 O sin 6 x (z) sin 6 2 (z) cos qp(z) 

- \ cos a Q cos 0 O cos 0 t (z) sin cp(z) - X sin a Q cos (z) sin# 2 (z) cos qp(z) 

+ \ sin a Q sin 6 ^ ( z ) sin qp( z) 

S 2 = \ cos a Q sin ^A 0 cos 6 2 (z)sincp(z) - cos a Q cos ^ 0 sin ^ (z) sin (z) sincp(z) 

+ cos a Q cos cos 0 l (z) cos cp(z) - sin a Q cos (z) sin# 2 (z)sincp(z) 

- sin a Q sin <9 X (z) cos cp(z) . 


Then the angle <A( z ) measured as previously defined is given by 

'A(z) = tan” 1 

With the derivation of the angles 0(z) and a(z) defining relative sun position, it is possible by 
four-point linear interpolation to determine the bending- moment components and the torsional- 
torque coefficient at all points along the boom length and hence solve the differential equation. 

It should be noted that if large-angle deflection is to be studied, the function <A(z) must be dif- 
ferentiated with respect to <p(z); this term must be introduced into the torsion equation as previ- 
ously discussed. 

METHOD OF SOLUTION 

Due to computer time limitations it is convenient to solve the thermal-bending-plus -twist 
problem in two parts. In part one, the numerical values of the elements of the matrices needed to 
define the temperature -dependent bending and torsion coefficients are determined. The output re- 
sults take a form that may be directly read into the computer as input data for part two. In part 
two the actual bending-plus -twist solution is performed and the desired results are computed 
outputs. 
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The following is a summary of the equations needed for solution of the problem and a descrip- 
tion of the particular technique used to generate the results of this report. 

The boundary value problem defining the temperature distribution, i.e., 

Kt d T / S ’ Z - = cr e T 4 (s, z) - J s a s T(s) [SL • n(s)] 
as 


for 0-si 2ttt 9 and 

Kt d2 ^ (S 2 ’ Z) = cr e [T 4 (s, z)-a s T 4 (s - 2ttt , z)] 

for 2ttt< s ^ p, where 


dT(s , z) 

ds |s-0 
s=p 


0 


is solved numerically by assuming a full set of initial conditions and matching-end conditions. This 
can be done with great accuracy on a digital computer by steepest-descent techniques. When the 
desired solution accuracy is achieved, the equation is solved again, together with the equations 
defining quantities that will be used in the bending-plus -torsion equations, that is, 


dT " (Z) - 1 T , , 

dT' - P T(s, z) 


dBM x (s, z) 
ds 


= e c EtT(s, z)y 2 (s) 


dBM y (s, z) 
ds 


= - e c Et T(s , z)x 2 (s) 


dQ(s, z) 
ds 


= p T(s, z) 


dV(s, z) 
<3s 


[e sin + - 2 ]t(s, z) - P^ r cos ( r 2 ) + Q( s ’ z ) 
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These quantities are evaluated at s = P and stored in a computer memory. The solution is repeated 
a finite number of times for values of 'P*, in the interval 


0 < V* < 2tt , 


and for values of a* in the interval 


. , rr 

0 < a* < ^ • 

A two-dimensional array of numbers is then obtained for each of the quantities 

T m (z), BM x (z), BM y (z) , and V( z ) . 

The function V(z) is further numerically differentiated with respect to V*; this, too, is stored in 
computer memory. 

The bending-plus-torsional equations are of the boundary value type; in order to determine a 
solution numerically by the Runge Kutta numerical integration techniques, a full set of initial con- 
ditions must be assumed. The boundary value problem is 


where 


d(9 i ( z ) 

dz 


BM X (z) 


BM y (z) 


cos 6. 


E I v cos cp(z) + gr sincp(z) 


dg 2 ( z ) 
dz 


BM X ( z ) BM y ( z ) 

" ET^“ sincp(z) + eT^“ cos <p(z) 


d 3 9(z) 


dz 3 


k 2 + \ 


‘ E e _c_f <?V(z) 
5cp(z) 


dtp 


p ( z ) 

dz 


e i (z) 


0 


< P ( Z )| 2 = 0 = 0 


6 2 ( Z )| 


z=0 


0 


dq)(z) 

' 3i_ z=0 


0 


and 


d 2 cp( z ) 
dz 2 


if tip is free to warp, 
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or 



0 , if tip is not free to warp. 


At each integration step the quantities (z) andi9 2 ( z ) are defined; hence, from the previously 
derived equations, 'A(z) and a( z ) are obtainable. With z ) and a( z ) given it is possible by four- 
point interpolation to determine the forcing functions BM X (z) , BM Y (z), and <?V(z)/dcp(z) for any 
relative sun position. 

The coordinates of any point on the boom in inertial space can be determined from the solution 
to the differential equations: 

dXj (z) 

37— - sin£ 2 (z) 


dY, (z) 

dz " - sin0, (z)cos 0 2 (z) 


dZ, (z) 

37- - cos 0 1 (z) cos e 2 (z) , 

where 

x . <•>!... = = Z .< ! >U = °- 

If large-angle deflection is anticipated, the torsional equation given above will no longer be 
valid, since the effect of bending on relative sun position has been assumed small. In order to cor- 
rect this, the following equation must be used in its place: 

d 3 <P(«) _ L + x tEe c r dV(z) ^(z) l dcp(z) 

dz 3 [ k K ^( z ) <M z )J ~ 

where d0(z)/^qp(z) is determined directly from the expression for </>(z) and dV(z)/d'P(z') is de- 
termined directly by four-point interpolation of the appropriate numbers in the matrix defining 
av(z)/a^*. 

SOLUTION OF BOUNDARY VALUE PROBLEM 

The simultaneous differential equations defined in the previous section are of the boundary 
value type. All boundary conditions except one are defined at the root z = 0. The one condition 
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defined at the tip is 


cp" (z) I = 0 , if tip is free to warp, 

I z — L 

cp ' (z) I = 0 , if tip is not free to warp. 

I 2 =L 

In order to obtain a solution for a boom of length L, the initial condition cp" (0) must be determined 
by a predictor-corrector or steepest- descent technique. It was found that this was not easily done; 
another approach had to be taken. 

The method used was simply to define an initial condition cp" (0) and carry the integration out 
with respect to z to some predefined point. The function cp" (z) (or cp' (z)) was then examined for 
zero crossing. For each value of z for which cp" (z) (or cp' (z)) was equal to zero, the boundary 
condition was satisfied and a solution was said to exist for that particular boom length. If for a 
particular value of cp" (0) , cp" (z) (or cp' (z)) was equal to zero at a number of points along z, the 
implication was that the same initial condition gave a solution for more than one boom length. 
Conversely, it is implied that if more than one boom length has the identical initial conditions, then 
a particular boom length can have a solution for more than one set of initial conditions. That is to 
say, for a boom of length L, more than one value of cp" (0) can satisfy the boundary condition at L. 

Thus, it can be concluded that even though the full set of initial conditions defines a particular 
solution uniquely, a particular boom of length L may have more than one deflected shape that will 
satisfy all boundary conditions. That is, for a particular sun orientation, a boom of given length L 
can have more than one static thermal equilibrium shape in inertial space. 


EXTENSIONS OF SOLUTION METHOD 

In deriving a means of solving the thermal -bending problem, a rather general method of attack 
for similar problems has been developed. In essence it is a method for determining the bending 
and twist of a thin-walled member of open section when the small-angle assumptions of Timoshenko T s 
problem are not valid over the entire member length relative to a fixed coordinate system. 

This method is particularly suited for the solution of problems where: 

1. the deflection of any short element of length is small-angle but the total deflection is large- 
angle relative to a fixed reference frame 

2. the bending is non-planar 

3. the position and orientation of every cross section relative to a fixed reference frame must 
be ascertained 

4. the forcing function is non-uniform along the length and may be position- and 
orientation-dependent 

5. the transverse-torsional coupling is expected to be a significant effect. 
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To apply this method we must consider the constraints of the problem and modify correctly 
the appropriate steps in the method. 

In most bending problems, the forcing function is usually either a definable stress distribution 
or a force distribution. 

If a stress distribution is the forcing function, equations describing the bending-moment com- 
ponents and torsional torque must be derived relative to the local reference frame in a manner 
similar to that presented. The equations of bending, twist, and transformation will be unchanged. 

If a force distribution is the forcing function, major modifications must be incorporated since 
the bending equations are no longer applicable as given. 

Assume that the forcing function can be described as 

F (z) = F x ( z ) i + F y ( z ) i • 


Then 


and 


d2 P x (z) 

dz 2 


F , <o 

E I 


d 2 Py (z) _ F y («) 

dz 2 F 

are the equation of bending. The derivative of the curvature components must be determined 
analytically, and the equations for the third derivatives of 0 1 (z) and 0 2 (z) solved. The torsion 
equation remains unchanged; however, the applied torque is given by the expression 

T sc CO = ( e + S ) x F(z) , 

where (e +?) is the vector between the shear center and the centroid of the cross section at z. Since 
the order of the system has been increased by four, four addition boundary conditions must be given. 
These are 


( 2 ) 


- P, 


(Z) l=L 


d£> (z) 


dz 


dPy (Z) 


z=L 


dz 


= 0 


if there is a zero tip constraint. 

It is obvious that the boundary value problem now is extremely difficult to solve; this must be 
done by a sophisticated steepest-descent or other technique. 
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FIGURES 7 THROUGH 13 APPLIED TORQUES, BENDING 
AND TORSION, VARIABLE OVERLAP 


A series of figures that show how the temperature varies about the perimeter of the boom will 
not tell much about the resultant bending and torsional effects of thermal stresses. Hence, figures 
of this type will not be provided. One question of interest, however, will be answered: what is the 
effect of a change in overlap angle < f>? Since boom overlap angle has a significant effect upon the 
temperature distribution, and hence the thermal bending, a series of curves are provided that show 
how the bending and torsional components vary with sun position for various overlap angles. 

Figure 7 shows the boom cross section, the direction of the X 2 (0), Y 2 (0) body axes, and the 
various labeled sun positions <A 0 to be studied. These sun positions are labeled as shown on the 
succeeding figures giving the bending and torsional components. These positions also correspond 
to the initial sun orientations for the various static -deflection solutions shown. 


Furthermore, only the case of counterclockwise seam orientation will be studied, since the 
coefficients and results associated with the other case differ only by \. 


Figure 8 shows the variation in the compo- 
nent of the bending moment about the X 2 (z) 
body axis as the sun position *P 0 is changed from 
0° to 360° for various boom overlaps (0°, 45°, 
90°, 135°, 180°). 
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Figure 8— Bending moment component about X 2 (z) 
body axis for boom overlap <$> vs. sun position i// Q . 



Figure 7 — Case I initial sun orientations to be studied. 


Figure 9— Bending moment component about Y 2 (z) 
body axis for boom overlap <£v$. sun position <// Q . 
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Figure 9 shows the variation in the component of the bending moment about the Y 2 (z) body axis 
as the sun position <P 0 is changed from 0° to 360° for various boom overlaps (0°, 45°, 90°, 135°, 180°). 

Figure 10 shows the change in magnitude and direction of the bending moment vector as <A 0 is 
changed from 0° to- 360° for various boom overlaps (0°, 45°, 90°, 135°, 180°). The five curves shown 
are obtained by plotting BM X (</<„) vs. BM Y (<£ 0 ) . The various sun positions labeled correspond to 
30-degree incremental changes in <p 0 and are also labeled on the curves. 


BENDING MOMENT 
Y- COMPONENT (lb. in. ) 




in. ) 




Figure 10— Variation in 


bending moment for boom overlap as sun positions change. 
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Figure 11 shows the angular amount of out-of -plane bending which would exist if the twist was 
assumed to be zero. For the zero-twist case, the temperature distribution is independent of boom 
length; the resultant bending moment due to the non-uniform thermal distribution would bend the 
non-symmetric beam in a direction perpendicular to its neutral axis. If the beam were to bend "in 
plane" it would bend away from the sun. Let <f(0 o ) be the angle measured in a right-handed sense 
from the positive X 2 (0) axis to the neutral axis for the case when the sun is at 0 O . By definition 
of the neutral axis. 


0A o ) 


i x bm y (-/-„) 

l y m x^o) ’ 


Let ¥ be the angle measured from the negative y 2 (0) axis in the direction of increasing s to the 
axis of "in plane" bending. If the sun is at ^ 0 , 


V = + 77 , 


where 


0 < V < 277 . 

Let 0 be the angle measured from the negative Y 2 (0) axis in the direction of increasing s to the 
axis of actual zero-twist bending (axis normal to neutral axis). It follows that 


0 



where 0 <0 < 2tt. The angular amount of out-of-plane bending for zero twist is shown by plotting, 
in Figure 11, V -0 vs. ^ 0 . The five curves shown correspond to the five different overlap angles 
considered. 



Figure 11— Angular amount of out-of -plane bending 
for zero-twist case vs. sun position. 



Figure 12— Variation of V(i/f Q ) vs. 
sun position '/'o f° r boom overlap. 
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Figure 12 shows the variation in the func- 
tion v(p, y^) as the sun position y^ is changed 
from 0° to 360°. It will be recalled that the 
derivative of this function with respect to lon- 
gitudinal arc length is proportional to the 
thermal torque under "Thermal Torque." The 
five curves shown correspond to the five over- 
lap angles considered. 

Figure 13 shows the variation of the 
derivative of the function v(p, \p Q ) shown in 
Figure 12 with respect to as sun position 
is changed from 0° to 360°. 


COMMENTS ON FIGURES 7 THROUGH 13 

Figures 7 through 13 provide a means of deriving a quantitative estimate of the effects of a 
change in overlap angle on the thermal bending moment and the thermal torque coefficient. By 
making use of the data provided on Figures 8, 9, or 10 and the equations 

i BM 

y / \ — 1 y 2 

X i ( z ) “ 2 ET z 

y 

1 BM x 

Y i ( z ) ~ 2 E I z2 ’ 

X 

a first-order approximation of the deflection associated with the trivial (cp" (0) = 0) zero-twist 
solution to the boundary value problem can be obtained. 

The non-trivial solution, however, cannot be simplified; and the data provided on these figures 
yield little information as to the actual thermal equilibrium shapes. Only after examination of a 
particular case can we gain some insight as to how to predict them, see below. 


SOLUTION TO A PARTICULAR PROBLEM 

The following is a list of the magnitudes of the geometrical and physical constants used to 
obtain a solution to the derived equations. The magnitudes stated characterize a silver-plated 
DeHavilland boom of the type now being considered for use as a Gravity Gradient boom on the ATS 
and jRAE satellites. 

J s = Solar radiation intensity = 3.065 BTU/(hr in 2 ) 
a s = Absorbtivity = 0.13 
P = Perimeter of cross section = 2" 
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e = Emissivity = 0.035 

K = Thermal conductivity = 6.0 BTU/(in F° hr) 
t = Thickness of cross section = 0.002 M 
e c = Thermal expansion coefficient = 0.104 X 10~ 3 4 * 
T 0 = Absolute ambient temperature = 535 R° 

E = Young's Modulus = 19 X 10 6 * * lb/in 2 
a - Radiation coefficient = 0.121 X 10" 10 BTU/(hr in 2 ) 
4> = Boom overlap angle = 100 degrees 
G = Shear Modulus = 6 x 10 6 lb/in 2 . 


Before a detailed discussion of each figure showing the results of the computer 
initiated, a brief mention of the limitations and implications of these results will be 


solution is 
given. 


1. It is felt that the equations defining the temperature distribution represent a worst-case 
estimate. That is, the actual temperature gradients will be smaller than those used in this analysis 
because 


a. internal radiation has been neglected 

b. zero contact in overlap region is assumed 

c. diameter changes are assumed negligible 

d. the boom is assumed homogenous 

e. actual boom oscillation in space will have an averaging effect. 


2. The boom is assumed to be perfectly clamped at the root. This condition in actual practice 
is impossible to obtain. It will be shown that, for a perfectly clamped boom under a given length 
and deflected by the predefined thermal stresses, only the trivial solution to the torsional part of the 
boundary value problem exists. This implies that short booms will be bent but not twisted by 
thermal stresses. Any twist, observed experimentally on short booms, is indicative of the non- 
perfect clamping and non-uniform physical characteristics of the boom itself. 


3. This analysis cannot and is not intended to give the exact deflected shape of an actual boom. 

However, it can and will bring to light many properties of thermal bending not generally known. It 

will be shown that some very unusual results can be obtained; for instance, part of the boom 

could actually corkscrew and bend back at the sun. Such cases are very special and involve 

simplifying assumptions; they are examples of what can happen rather than what will happen in 

flight. 
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In order to transmit the maximum information obtained from this analysis with the least con- 
fusion, the number of figures presented has been kept to a minimum. 

As was stated in the body of the report, the method used to obtain a solution was to pick arbitral 
values of the free unknown initial condition cp" (0), run the program, and observe at what values of 
arc length z the boundary condition is satisfied. These points are recorded, <p" (0) is changed 
slightly, and the analysis is repeated. At each solution point, four quantities are recorded: 

1. The initial condition 

2. The twist at this solution point, "the tip" 

3. The magnitude of the deflection at the tip 

4. The direction of the tip deflection (bending). 

The figures provided show the results of computer runs over a full range of initial conditions for 
a boom with a tip free to warp; that is, qp" (L) = 0. For the problem when the tip is not free to 
warp, qp' (L) = 0, the analytical solutions are not grossly different; in fact, the same general state- 
ments can be made about both cases. Hence the latter set of figures will not be included. 

The four quantities listed above are each plotted against boom length (solution length), for 
eight different values of sun position; this gives 32 figures: Figures 14a through 21d. Table 1 
shows the relation of figure number and letter to quantity plotted and sun position. 

Table 1 


Guide to Figures 14a Through 2 Id. 


't’o 

Sun Position 
(degrees) 

<P" (0). 
Initial 
Condition 

<P(L), 

Tip 

Twist 

R(L), 

Tip 

Deflection 

0(L). 

Direction 
of Tip 
Bending 

0 

14a 

14b 

14c 

14d 

50 

15a 

15b 

15c 

15d 

90 

16a 

16b 

16c 

16d 

130 

17a 

17b 

17c 

17d 

180 

18a 

18b 

18c 

18d 

230 

19a 

19b 

19c 

19d 

270 

20a 

20b 

20c 

20d 

310 

21a 

21b 

21c 

2 Id 



FIGURE NUMBERS 
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R(L), TIP DEFLECTION (inches) ¥>"(0), INITIAL CONDITION 




0 50 100 150 


Figure 14c— Tip deflection R(L) vs. L, SOLUTION LENGTH (feet) 

solution length L for \p 0 = 0. 

Figure 14d— Direction of tip bending 6 (L) vs. 
solution length L for <// Q = 0°. 
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R(L), TIP DEFLECTION ( inches ) V>"(0), INITIAL CONDITION 



Figure 15a— Initial condition cp" (0) vs. 
solution length L for = 50°. 




Figure 15b— Tip twist qp(L) vs. 
solution length L for \p Q = 50°. 



Figure 15d— Direction of tip bending 0(L) 
solution length L for \p Q = 50°. 


Figure 15c— Tip deflection R(L) vs. 
solution length L for = 50°. 







R(L), TIP DEFLECTION (inches) ^"(O), INITIAL CONDITION 


L, SOLUTION LENGTH (feet) 



'(2) (3P (4) (5) 


L, SOLUTION LENGTH (feet) 



(2) (3) (4) 


Figure 18a— Initial condition cp" (0) vs. 
solution length L for \ p Q = 180°. 


Figure 18b— Tip twist cp(L) vs. 
solution length L for \p Q = 180°. 



L, SOLUTION LENGTH (feet) 

Figure 18c— Tip deflection R(L) vs. 
solution length L for i// Q = 180°. 


PLANE L, SOLUTION LENGTH (feet) 

Figure 18d— Direction of tip bending 0(L) vs 
solution length L for t// n - 180°. 






R(L), TIP DEFLECTION (Inches) " (0), INITIAL CONDITION 


270 



Figure 20a— Initial condition cp" (0) vs. 
solution length L for = 270°. 


75 


(4) 



0 50 100 150 


L, SOLUTION LENGTH (feet) 

Figure 20c— Tip deflection R(L) vs. 
solution length L for t// Q = 270°. 



Figure 20b— Tip twist <p(L) vs. 
solution length L for 0 O = 270°. 



L, SOLUTION LENGTH (feet) 


Figure 20d— Direction of tip bending <9(L) vs. 
solution length L for — 270°. 
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R (L) TIP DEFLECTION (inches) INITIAL CONDITION 



Figure 21a— Initial condition q>" (0) vs. 
solution length L for \p Q = 310°. 



L, SOLUTION LENGTH (feet) 


Figure 21c— Tip deflection R(L) vs. 
solution length L for \p Q = 310°. 


Figure 2 1 d — Direction of tip bending 6?(L) vs. 
solution length L for 0 O = 310°. 



Figure 21b— Tip twist <p(L) vs. 
solution length L for^ 0 = 310°. 



50 100 

L, SOLUTION LENGTH (feet) 


150 
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RESULTS OF DIGITAL SOLUTION, FIGURES 14 THROUGH 25 


Since different results will be obtained for different sun positions, eight different positions 
are examined; these correspond to those positions labeled on Figure 7. Figures 14a through 14d 
correspond to sun position 1, Figures 15a through 15d correspond to sun position 2, and so on up 
to sun position 8. For each sun position four different quantities are plotted vs. boom length. 

Figures 14a, 15a, etc., plot initial condition q>" (0) vs. boom length L for such positions 1, 2, 
etc., respectively. Each point of the curves plotted on these figures corresponds to a point at which 
the boundary condition is satisfied for the assumed initial conditions q>" (0) . Hence, if a particular 
boom of given length is to be analyzed, the initial condition or conditions which will yield a static- 
thermal -equilibrium solution for that boom and sun position can be determined from the appropriate 
figure. From these curves it is apparent that for most boom lengths two of the initial conditions 
that yield a solution are easily definable; however, all of the other initial conditions that also yield 
solutions lie in the asympotic region of the curves and are extremely difficult to obtain. 

Figures 14b, 15b, etc., plot the magnitude of the tip twist cp(L) vs. boom length L for sun posi- 
tions 1, 2, etc., respectively for each solution shown in the M a M series of figures. Each point of 
the curves plotted gives the magnitude of the twist at the boom tip but says nothing about what the 
twist will be between the root and the tip. Hence, for a particular boom being analyzed it is pos- 
sible to determine what the magnitude of the tip twist will be for each of the booms’ thermal- 
equilibrium shapes. 

Figures 14c, 15c, etc., plot the magnitude of the tip deflection R(L) vs. boom length L for sun 
positions 1, 2, etc., respectively, for each solution shown in the ”a” series of figures. The magni- 
tude R(L) is measured in inertial space and is determined from the equation 


R(L) = [x, (L ) 2 + Yj (L) 2 ] 1/2 , 

where X x (L) andYj (L) are the components of the boom tip in the (x 1? Y x ) plane, fixed in inertial 
space. Hence, for a particular boom being analyzed it is possible to determine the tip deflection 
for each of the booms’ thermal -equilibrium shapes. 

Figures 14d, 15d, etc., plot the direction of the tip deflection 9(h) vs. boom length L for sun 
positions 1, 2, etc., respectively, for each solution shown in the "a” series of figures. The magni- 
tude of 9(h) is 


9(h) 


77 

"2 + 


tan ’ 1 


(L) 

*i (L) * 


Figure 22 shows the geometry of these relations. In each of the ”d” series of figures, two 
marks are shown on the vertical axis. These marks correspond to the direction of bending if the 
boom were to bend ”in plane” and the direction of bending if the boom were to bend with ’’zero 
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Figure 22— Geometry of tip deflection. 


twist." Hence, for a particular boom being analyzed it 
is possible to determine the direction of the tip deflection 
for each of the boom T s thermal -equilibrium shapes. It 
should be noted, however, that these figures give no infor- 
mation as to how the direction varies along the length of 
the boom. 

In each set of four figures corresponding to a partic- 
ular sun orientation the curves appearing in each figure 
are numbered. Points adjacent to the same number on the 
four figures correspond to the solution generated from the 
same set of boundary conditions. 


From a close examination of the preceding Figures 14a through 21d, certain general facts 
about thermal bending plus twist become immediately apparent. 

1. For a particular sun position, a boom of a given length may have more than one thermal- 
equilibrium shape. 

2. The different thermal-equilibrium shapes correspond to the different initial conditions that 
satisfy the boundary condition of the boom under study. 

3. For a particular sun position and given boom length, the magnitude and direction of the tip 
deflection can be radically different for different thermal- equilibrium shapes. 

4. The means to determine the region of stability for each of the thermal- equilibrium shapes 
is not apparent from this analysis. Hence, the most probable equilibrium shape cannot be 
determined. 

5. The inclusion of transverse-torsional coupling through the cross-sectional-orientation de- 
pendence of the thermal-stress distribution in the six-dimensional analysis yields results 
that cannot be predicted from either a two- or a three-dimensional analysis. 

6. If the temperature-distribution equation could be replaced by a more accurate approxima- 
tion results, #1, 2, 4 and 5 would remain unchanged; however, if the thermal gradients are 
significantly reduced, result #3 would most probably have to be tempered. 


In order to emphasize the fact that a given boom may have more than one thermal-equilibrium 
shape for a given sun position, three additional sets of figures are included. On each of the figures, 
curves are shown that correspond to some of the thermal -equilibrium shapes of a one-hundred foot 
boom. The following is a list of the sun positions t£ 0 and initial conditions <p" (0) used to compute 
these curves: 

'K <p" (Q) 

Figure 23 

Curve (1) 90° 0.107526 X 10' 3 

Curve (2) 90° 0.31657 x 10' 4 
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Figure 24 

'P 

l_o 

<p"(0) 

Curve (1) 

130° 

0.77723 X 10" 3 

Curve (2) 

130° 

-0.23172 X 10~ 4 

Curve (3) 

130° 

-0.80193 X 10' 4 

Curve (4) 

130° 

-0.8034491 X 10~ 4 

Figure 25 

Curve (1) 

230° 

-0.14515 X 10' 5 

Curve (2) 

230° 

-0.3824 X 10‘ 4 

Curve (3) 

230° 

-0.125083 X 10 -3 . 


By including only these few thermal- equilibrium shapes it is not implied that these are the only 
shapes or the most probable shapes, but that they are the shapes for which the appropriate initial 
conditions can be determined from the preceding figures. The other equilibrium shapes have 
initial conditions lying in the asympotic regions and are extremely difficult to determine. 

The sun positions studied in these figures were chosen because the resultant thermal- 
equilibrium shapes vividly illustrate that they are not even approximately equal and hence cannot 
be approximated by any simple function. 

Each set of Figures 23, 24, and 25 has four individual figures associated with it labeled a, b, 
c, and d. Table 2 shows the relation of figure number and letter to subject matter and sun position. 

The ,T a TT figure of each set shows the projection of the boom T s thermal-equilibrium shape on 
the (X x , Y t ) inertial plane. That is, the coordinates X 1 (z) vs. Y t (z) are plotted for each of the 
derived thermal-equilibrium shapes. The shapes shown are not the only ones that exist but they 
are the most easily definable. 

The "b TT figure of each set shows how the twist cp(z) varies along the boom length. The curves 
labeled 1, 2, etc., correspond to those similarly labeled in the a, c, and d figures. 


Table 2 

Guide to Figures 23a Through 25d. 


'/'o- 

Sun Position 
(degrees) 

Projection of 
Thermal- Equilibrium 
Shape 

Twist vs. 
Arc Length 

Tip 

Deflection vs. 
Arc Length 

Bending 
Direction vs. 
Arc Length 

90 

23a 

23b 

23c 

23d 

130 

24a ■ 

24b 

24c 

24d 

230 

25a 

25b 

25c 

25d 



FIGURE NUMBERS 
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Yj ( inches) 


URVE <P"(0) 

1 -0.77723 X 10' 3 

2 -0.23172 X 10‘ 4 

3 - 0.80193 X 10' 4 

4 - 0.803449 X 10' 4 



ction of thermal-equilibrium 
inertial plane for L = 100 ft. 


( 2 ) 


8(z), BENDING DIRECTION (degrees) -r, *(*)/ TWIST (degrees) 





TIP DEFLECTION (inches) 



Figure 25a — Projection of thermal-equilibrium 
shape in (X v Yj) inertial plane for L = 100 ft, 
= 230 ° * 



Figure 25c— Tip deflection R(z) vs. arc length z for 
thermal-equilibrium shapes shown in Figure 25a. 



Figure 25b— Twist cp(z) vs. arc length z for thermal- 
equilibrium shapes shown in Figure 25a. 



z, ARC LENGTH (feet) 


Figure 25d— Bending direction 0 ( z) vs. arc length ; 
for thermal -equilibrium shapes shown in Figure 25a 



The TT c” figure of each set shows how the magnitude of the deflection increases along the boom 
length. The magnitude R(z) is derived from the equation 


R ( z ) “ [ x i ( z ) 2 + Y i ( z ) 2 ] 1/2 • 

The "d" figure of each set shows how the direction of bending changes along the boom length. 
The quantity 0 { z) which is plotted is defined by the equation 


<9(z) 



where 0 < 0(z) < 2tt. 

Examination of Figures 23, 24, and 25 discloses certain facts about thermal bending: 

1. For this particular boom and temperature distribution, the thermal -equilibrium shapes are 
grossly different from each other. 

2. Each thermal-equilibrium shape corresponds to a distinct pattern of twist. 

3. The thermal-equilibrium shapes that correspond to twist patterns of low magnitude are 
nearly planar in bending. 

4. For the cases studied it is apparent that the direction of bending has its most significant 
change in the interval of length along which the twist is going through its first half cycle. 

5. The large out-of-plane bending, evident for some thermal- equilibrium shapes, is a result 
of the coupling between bending and torsion. This coupling over the length is a result of the 
fact that the direction of bending at any point depends not only on the direction of the 
thermal-bending-moment vector at the point in question but also on the shape of the boom 
up to this point. 

6. The direction of deflection depends on more than just the local stress distribution and 
position. This can be seen from an analysis of the equations of bending. Since the order 
of this system of equations is eight rather than six, it follows that the six coordinate mag- 
nitudes at a point and the forcing function are not enough to define the coordinate magni- 
tudes at the point an infinitesimal distance away. To determine these magnitudes, two 
additional conditions must be given which characterize the boom's shake between the root 
and the point in question. 

7. The deflected shape of the boom up to a particular point affects the deflected shape after 
that point, and the thermal stresses can induce a significant amount of out-of -plane bending. 
These two effects can be combined through the transverse torsional coupling. This may 
bring about much more out- of- plane bending than would be expected from considering 

only the resultant bending- moment variation with sun position as shown in Figure 10 or 
the results shown in Figure 11. 
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ANALYTIC SOLUTION TO TORSION EQUATION 


The torsion equation as previously derived; i.e., 


fc t Ee c r dW(z) 3Wz)l , 

( z ) [C t + *- Cj ^(z) <9<p(z ) J'P ( 2 > " 0 


is highly nonlinear and coupled with the bending equation for large-angle twist. If, however, one 
restricts the discussion to small-angle twist and bending, 

*(*) = 'Z', - \q>(z), 


and the torsion equation can be written in the form 


qp'" (z) - [A 2 + Bqp(z) + Dcp 2 (z)]qp'(z) - 0, 

where A 2 , B, and D are constants whose magnitude can be determined from the coefficient of the 
derived torsion equation given above and from a parabolic representation of aV(z/^W(z) about the 
initial sun orientation (see Figure 13). 

When the above equation is integrated once, it can be put in the form of a general elliptic 
equation; that is, 


qp " (z) - cp" (0) + A 2 qp( z ) + 3 qp 2 (z) + 3 qp 3 (z) , 

where cp" ( 0 ) is the constant of integration equal to cp" ( z ) I when qp(z)| =0. 

z = 0 } z= 0 

In Reference 5 it is shown that by an additional integration this equation can be put in the form 

cp' (z ) 2 = cp' (0 ) 2 + 2qp>" (O)cp(z) + A 2 qp 2 (z) + 3 cp 3 (z) +3 cp 4 (z) 

and that a transformation 


£ = £W 


can be found which further reduces the equation to the form 


(S) ! = • 
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This equation is well known and the solution to it is given by 


€ “ sn(z, m) , 

where sn(z, m) is the elliptic sine of Jacobi. Inverting the transformation <f = <f(<p) gives the equa- 
tion defining cp(z); it is given in terms of the Jacobi elliptic functions. 

Function cp(z) can be computed for any particular case; but a general solution cannot be obtained, 
since the solution requires the determination of the roots of a quartic equation and the magnitude of 
<p" (0). The quantity 9 " (0) is a function of the boundary condition and must be obtained from the 
solution of a transcendental equation written in terms of the Jacobi elliptic functions. 

Since a general solution to the elliptic equation cannot be found, the more restrictive assump- 
tion that 


where 


A 2 » Bcp(z) + Dtp 2 (z) 


A 2 


C-tEe 

C 


av(^) 

JH'(z) 


must be made. It follows from the parametric magnitudes listed on page 36 and the equations given 
on page 16 that 


C = 0.032 lb in 2 

Cj = 1413.7 lb in 4 

t Ee r = 0.0985 lb , 

C 


and hence 


A 2 = 


1413.7 


0.032-0.0985 


*V(^o) 

d^z) 


The actual magnitude of dW^P 0 )/d^(z) can be obtained approximately from the curves provided in 
Figure 13. From the shape of these curves it is apparent that the above assumption is reasonable 
in the interval 


50° < % < 200° . 

Outside of this interval the assumption is valid only for very small <p(z). 
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In order to interpret the numerical solutions shown in Figures 14 through 21 and extrapolate 
these results to similar problems, it is useful to obtain the general solution of the linearized torsion 
equation for both positive and negative A 2 and for perfect and imperfect clamping. 

It follows that the solution to the linearized torsion equation 


cp'" (z)~A 2 cp'(z) - 0, 


subject to the boundary conditions 


<P(0) ~ <p" (L) - 0 , 


is: 


Case I 


A 2 > 0 , <p' (0) / 0 , 

cp ' (0) f sinhA(L-z)\ 

<P( Z ) “ A ~ ^tanh AL Eosh AL /’ 
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Case II 


A 2 > 0 , cp' (0) = 0 


<p(z) - 0 , 


Case III 


A 2 < 0 , <p' (0) / 0 , 


cp(z) 


V' ( 0 ) 


{ 


sin A(L - z) 
tanAL cos AL 


} 


Case IV 


A 2 < 0 , cp' (0) = 0 , L = 2A ( 2n “ !) 


cp " ( 0 ) 

cp(z) ^ (1-cosAz) , 


Case V 


A 2 < 0 , cp' (0) = 0 , L / 2 H ( 2n " l ) 


where 


<p(z) - 0 , 


A = A 2 


1/2 


ii i 


i i 


i i 



A comparison of these equations with the solution shown in Figures 14 through 21 yields a 
number of results which further explain the observed phenomena. These can be summarized as 
follows: 

1. For the particular case of perfect clamping, q>' (0) = 0, a non-trivial linear small-angle 
solution exists only if A 2 is negative and the boom length L is given by 

L = -jf (2n-l) . 


In the sun interval 


Figure 13 shows that 


50° < ^ < 200° . 


* V (*0> ■ , , . 2 

d*V(z) 13 in 


and hence 


A 


0.032 - 0.0985 x 1. 5 
1413.7 


1/2 


= 0.009049 in' 1 , 


or, for n = 1, 


L = 3S - 173.6" 


14.5' . 


Figures 16 through 18 correspond to sun orientations W 0 within the interval 


60° < ^ < 200° . 

In these figures the numerical solution shows that only booms of specific lengths will have non- 
trivial small- angle solutions and that these solution lengths can be approximated by the equation 


L = (2n - 1) x 14.5 ft n = 1, 2, - • ■ 


as predicted by the analytic solution Case IV and Case V. 
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Figures 15 and 19 correspond to sun orientations v^ 0 of 50 and 230 degrees, respectively. For 
both cases A 2 = 0; hence, the linearizing assumptions made are not valid. The numerical solution, 
however, shows that the elliptic equation can have a small-angle solution for any length boom 
greater than 14.5 ft. 

Figures 14, 20, and 21 correspond to sun orientations of 0, 270, and 310 degrees, respec- 
tively. For these orientations A 2 is positive and the linearizing assumption holds for very small- 
angle twist. The numerical solutions shown bear out the analytically predicted result (Case n) that 
a non-trivial small angle solution does not exist. 

2. For the case of imperfectly clamped booms of any length, if the twist is such that the 
linearizing assumptions are valid over the entire length, a non-trivial small angle twist solution 
can be obtained. Then the magnitude of the twist is proportional to the magnitude of the clamping 
imperfection. 

3. From the equation defining A 2 ; i.e.. 


A 2 


_l 

c , 


C-t Ee 

C 




it is apparent that the relative magnitude of the torsional rigidity and the thermal torque coefficient 
dictates whether the twist will be hyperbolic or trigometric over any interval for which the linear- 
izing assumptions hold. This equation implies that if the torsional rigidity C can be significantly 
increased over the stated value herein, A 2 will be a positive constant approximately independent of 
sun orientation; the twist will therefore be hyperbolic over each interval. 

Furthermore, if A 2 is positive and approximately independent of sun orientation, the torsional 
solution will be 


<p(z) 



*|tanh AL - 


sinh A(L - z ) ^ 
cosh AL J 


and the small-angle restriction no longer need be made. This implies that, for the case of 
C»tEe c r > * or ^ thermal equilibrium shape is unique and the twist is directly 

proportional to the clamping imperfection. 

It is interesting to compare the very different torsional rigidities of a seamless cylinder and 
a cylinder of open section, both having the same cross-sectional properties. 

For the case of the seamless cylinder, the torsional rigidity is given by 

C = t Pr 2 G 
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and, as previously stated, the torsional rigidity of the cylinder of open section is 

C = J t 3 PG . 

Thus, the torsional rigidity of a seamless cylinder is 3r 2 /t 2 times greater than that of a cylinder 
of open section, which for the constants of this problem is equal to 0.3675 x 10 s . 

Hence, it is not unreasonable to expect that a Gravity Gradient boom could be designed that 
would have an effective torsional rigidity significantly larger than the coefficient of the thermal 
torque and hence have a unique thermal-equilibrium shape for every sun orientation. 


CONCLUSIONS 

The most important facts shown in this analysis of a particular type of Gravity Gradient boom 
may be summarized as follows: 

1. A method has been developed by which the non-planar thermal bending plus twist of very 
long thin-walled members of open section may be studied. 

2. The thermal-equilibrium shape of a Gravity Gradient boom in a solar thermal field is not 
unique. 

3. The various possible thermal- equilibrium shapes may radically differ from each other in 
magnitude, direction, and shape. 

4. Transverse -torsional coupling is a very significant effect in booms having high thermal 
gradients and low torsional rigidity. 

5. A means of performing a stability analysis on the various thermal -equilibrium shapes is 
not apparent from this analysis. 

6. A unique function that can be used in the dynamic analysis of a Gravity Gradient satellite to 
approximate the thermal bending of the booms does not exist for the particular case studied. 

7. Thermally induced twist cannot be assumed to have small magnitude. 

8. The thermal- equilibrium shapes that have twist of low magnitudes associated with them 
tend to bend in a nearly predictable direction. 

9. For every sun orientation studied herein the torsional rigidity of 0.032 lb in 2 is large 
enough to prevent booms of length less than 14.5 ft from assuming more than one thermal- 
equilibrium shape. For longer length booms, however, more than one thermal- equilibrium 
shape can satisfy the boundary conditions of the problem. 

10. The direction of deflection of any point along the booms* length depends not only on the 
local conditions (that is, coordinate magnitude and stress distribution) but also on the de- 
flected shape between the root and the point in question. 
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11. It follows from conclusion 10 that the maximum amount of out-of-plane bending cannot be 
predicted from a simple study of Figure 11 which plots zero-twist out-of-plane bending vs. 
sun position. 

12. If the torsional rigidity can be significantly increased above the thermal torque coefficient, 
the number of thermal-equilibrium shapes will be reduced to one unique thermal- equilibrium 
shape. 

ACKNOWLEDGMENTS 

The author gratefully acknowledges the assistance of Mr. Benjamin Zimmerman, Head of the 
Analytical Studies Section, for the support he provided throughout the development of this report 
and Miss Kathryn Dohm for the typing of the original manuscript and its succeeding revisions. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland, April 5, 1967 
630-12-02-01-51 


REFERENCES 


1. Timoshenko, S., "Strength of Materials Part n," Princeton, N. J. : Van Nostrand, 1956. 

2. Timoshenko, S., "Theory of Bending, Torsion, and Buckling of Thin Walled Members of Open 
Section," Journal of the Franklin Institute , 239(3, 4, 5), 1945. 

3. Boley, B. A., and Weiner, J. H., "Theory of Thermal Stresses," New York: John Wiley, 1960. 

4. Southwell, R. V., "An Introduction to the Theory of Elasticity for Engineers and Physicists," 
London: Oxford University Press, 1941. 

5. Davis, H. T., "Introduction to Nonlinear Differential and Integral Equations," New York: 

Dover Publications, 1962. 


60 


N Ab A - Langley , 1967 


19 



fr The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 
important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless of 
importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distribu- 
tion because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Scientific and technical information generated 

under a NASA contract or grant and considered an important contribution to 
existing knowledge. 

TECHNICAL TRANSLATIONS: Information published in a foreign 
language considered to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information derived from or of value to NASA 

activities. Publications include conference proceedings, monographs, data 
compilations, handbooks, sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION PUBLICATIONS: Information on tech- 
nology used by NASA that may be of particular interest in commercial and other 
non-aerospace applications. Publications include Tech Briefs, Technology 
Utilization Reports and Notes, and Technology Surveys. 


Details on the availability of these publications may be obtained from: 


SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



